Dijet production as a centrality trigger for pp collisions at CERN LHC 
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We demonstrate that a trigger on hard dijet production at small rapidities allows to establish a 
quantitative distinction between central and peripheral collisions in pp and pp collisions at Tevatron 
and LHC energies. Such a trigger strongly reduces the effective impact parameters as compared to 
minimum bias events. This happens because the transverse spatial distribution of hard partons (x > 
f 0~ 2 ) in the proton is considerably narrower than that of soft partons, whose collisions dominate the 
total cross section. In the central collisions selected by the trigger, most of the partons with x > 10 -2 
interact with a gluon field whose strength rapidly increases with energy. At LHC (and to some extent 
already at Tevatron) energies the strength of this interaction approaches the unitarity ("black- 
body") limit. This leads to specific modifications of the final state, such as a higher probability of 
multijet events at small rapidities, a strong increase of the transverse momenta and depletion of the 
longitudinal momenta at large rapidities, and the appearance of long-range correlations in rapidity 
between the forward/backward fragmentation regions. The same pattern is expected for events with 
production of new heavy particles (Higgs, SUSY). Studies of these phenomena would be feasible 
with the CMS-TOTEM detector setup, and would have considerable impact on the exploration of 
the physics of strong gluon fields in QCD, as well as the search for new particles at LHC. 

PACS numbers: 12.38.-t, 13.85.-t, 14.80.Bn, 25.75.Nq 



I. INTRODUCTION 

The differences between peripheral and central colli- 
sions play a crucial role in the physics of heavy ion colli- 
sions. In pp collisions, a similar distinction can be made 
at the energies of the LHC and the Tevatron. This is pos- 
sible because of the appearance of two separate transverse 
distance scales at high energies. On one hand, as pre- 
dicted by Gribov Q, the essential impact parameters in 
hadron-hadron collisions increase with the energy. This 
has been observed e.g. in numerous experiments in elas- 
tic pp scattering, see Ref. for a review. On the other 
hand, the transverse spatial distributions of hard partons 
(with finite light-cone fraction, x) in the colliding nuclc- 
ons is only a weak function of x. For the gluon distri- 
bution this has been verified experimentally in studies of 
the t-dependence of photoproduction of heavy quarkonia 
off the nucleon jj] ■ The two scales allow us to classify pp 
collisions at collider energies. A schematic illustration of 
this idea is given in Fig.^ In collisions with large impact 
parameters there will be essentially no overlap between 
the hard partons (Fig.^). Only partons with x <C 1CT 2 
will overlap with significant probability. These peripheral 
events constitute a significant (in fact, dominant) part of 
the total inelastic cross section. The production of high 
p± jets, however, as well as of heavy particles, will be 
strongly suppressed. At small impact parameters, how- 
ever, the distributions of hard partons in the two colliding 
nucleons will overlap, and the probability of hard interac- 
tions will be greatly enhanced (Fig. QJ)) . This difference 
between the physics of soft and hard QCD processes (i.e., 
with x±,X2 of the colliding partons > 10 -2 ) gives us an 
opportunity to distinguish quantitatively between central 
and peripheral collisions at collider energies. 



Specifically, we propose here to use the production of 
(one or more) hard dijets near zero rapidity as a "central- 
ity trigger" for pp and pp collisions at Tevatron and LHC 
energies. At the LHC, such a trigger could be imple- 
mented with any of the central detectors. It will lead to 
a significant enhancement of the production of hadrons 
at small rapidities and drastic changes of the pattern of 
forward production, which could be probed, for example, 
by the TOTEM detector in combination with the CMS 
detector. Since the production of heavy particles, such 
as the Higgs boson or supersymmetric particles, is also 
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FIG. I: Schematic illustration of the two classes of pp colli- 
sions at high energies. The transverse spatial distributions of 
the hard partons (a; > I0 -2 ) is indicated by the dark shaded 
disks, those of the soft partons (x <C 10 -2 ) by the light shaded 
disks; 6 denotes the impact parameter of the pp collision, (a) 
At large b no overlap between hard partons occurs, (b) At 
small b the distributions of hard partons overlap, leading to 
production of hard dijets (and, possibly, heavy particles). 



greatly enhanced for central collisions, such a program 
could have considerable impact on the searches for new 
particles at the LHC. 

Another important application of the proposed "cen- 
trality trigger" is in the investigation of the physics of 
strong gluon fields in QCD. Numerous measurements of 
small-x phenomena at HERA have confirmed the fast 
increase of the gluon density in the proton predicted 
by perturbative QCD; for a review and references see 
Ref. 0. At LHC energies, and to some extent already at 
the Tevatron, the gluon density becomes so large that the 
interaction of high p± partons with the gluon "medium" 
becomes strong and multiple scattering effects cannot be 
neglected, for a review and references see Ref. Q- The 
new phenomenon one encounters here can be described 
as the breakdown of the Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi (DGLAP) approximation caused by the 
approach to the unitarity limit(Black Body Limit, BBL), 
in which partons with virtualities below a certain value 
interact with the other proton with the maximal strength 
allowed by s-channel unitarity. More quantitatively, a 
parton in one proton with longitudinal momentum frac- 
tion xji and virtuality p\ resolves partons in the other 
proton with 
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where s is the invariant energy squared of the pp colli- 
sion, see Fig. [21 In particular, l&rge-XR partons resolve 
small-x partons in the other proton. At the LHC en- 
ergy, y/a = 14000 GeV, for p x ~ 2 GcV/c and x R ~ 10~ 2 
one obtains x ~ 10~ 5 . Under this condition the inter- 
actions with gluons can approach the BBL, see Ref. p| 
and references therein. Present data on heavy quarko- 
nium photoproduction indicate that the gluon density at 
small x is maximum in the transverse center of the proton 
• Since the distribution of large-Xi? partons is likewise 
concentrated at small transverse distances, it is evident 
that the chances of approaching the BBL are maximum 
for central pp collisions, which can be selected with the 
proposed "centrality trigger" . 

The approach to the BBL in central pp collisions leads 
to a significant change of the initial and final state inter- 
actions in the hard QCD processes. The interactions with 
the "high density gluon medium" suppress the spectrum 
of low transverse momentum partons and enhance the 
high momentum tail. Typically, all partons which are not 
involved in the hard interactions at the scale much higher 
than the BBL scale, jj^.bbl, receive transverse momenta 
~ P±.bbl and experience a significant momentum loss 
due to the gluon radiation (a schematic illustration is 
provided in Fig. A forward would-be spectator com- 
ponent of the proton wave function is "pulverized" com- 
pletely, loosing its coherence. As a result the number of 
particles with xf > 0.1 in the proton fragmentation re- 
gion will strongly diminish, while the average transverse 
momenta of these particles will grow to values compa- 
rable to p_l,bbl- Much more energy will be released at 
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FIG. 2: Schematic illustration of the effect of the black- 
body (unitarity) limit on hardon production in the for- 
ward/backward rapidity region in central pp collisions. A 
small-a; spectator parton (i.e., a parton not involved in the 
centrality trigger) from the left proton propagates through the 
strong gluon field (indicated by the shaded area), acquiring a 
large transverse momentum, p_l,bbl 3> Aqcd- The small-a; 
parton is then resolved in a collision with a large-x_R parton 
from the right proton, resulting in hadron production in the 
backward rapidity region. 



small rapidities as compared to minimal bias events. 

Although the phenomena discussed here are higher- 
twist corrections to the inclusive cross section and the 
transverse spectra of the production of sufficiently heavy 
particles (Higgs, SUSY) , in the vicinity of the BBL they 
will strongly modify the overall structure of the final 
states. In particular, they will change the pattern of 
radiation of moderate p± jets and the Sudakov form fac- 
tors for dijet production. Thus, the understanding of 
these phenomena is important also for effective searches 
for new physics at the LHC. 

The basic idea of the proposed "centrality trigger" is 
that the restriction to events with production of a hard 
dijet strongly reduces the effective impact parameters in 
high-energy pp and pp collisions as compared to mini- 
mum bias events. One can further narrow the distribu- 
tion of impact parameters by requiring the presence of 
multiple dijets in the same event (this was first studied 
within the framework of a multiple interaction Monte 
Carlo model in Ref. 0). The actual reduction which 
can be achieved in this way depends, to some extent, on 
possible spatial correlations of hard partons in the trans- 
verse plane. In this respect we make a surprising obser- 
vation, namely that the data on double dijet production 
in pp collisions at the Tevatron obtained by the CDF Col- 
laboration indicate significant spatial correlations of 
partons in the transverse plane. 

The paper is organized as following. In Section ITT1 we 
review the available information on the impact parameter 
distribution in the generic inelastic pp collisions. In Sec- 
tion IIIII we summarize our knowledge of the transverse 
spatial distribution of gluons in the nucleon and study 
its dependence on the resolution scale. In Section llVl we 
calculate the impact parameter distribution in pp colli- 
sions with production of a hard dijet near zero rapidity, 
and compare it with the impact parameter distribution 
for generic inelastic collisions. We demonstrate that hard 
dijet production acts as a "centrality trigger". We also 
discuss the extension to production of multiple dijets and 
the role of possible correlations in the transverse spatial 
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distribution of gluons. In Section [V] we investigate the 
role of the "centrality trigger" in approaching the BBL 
in central pp collisions. We show that the trigger es- 
sentially eliminates collisions at large impact parameters 
where the soft interactions are not black. In Section IVT1 
we list the novel characteristics of the final state in central 
pp collisions which follow from the proximity to the BBL 
of the spectator parton interactions at small transverse 
distances. Our conclusions are presented in Section lYlII 



II. IMPACT PARAMETER DISTRIBUTION 
FOR GENERIC INELASTIC pp COLLISIONS 

We begin by summarizing the available information 
about the impact parameter distribution of the cross sec- 
tion for generic inelastic pp collisions at high energies. 
Most of our knowledge here comes from pp elastic scat- 
tering, which has been studied in numerous experiments 
see Ref. for a review. By unitarity (i.e., the optical 
theorem) the pp elastic amplitude contains information 
also about the total (elastic plus inelastic) cross section, 
and thus about the inelastic cross section. 

It is well known from studies of pp elastic scattering 
that the radius of strong interactions (the average impact 
parameter) increases with the collision energy, s. The t— 
slope of the elastic cross section, B, grows as 



B(s) = B(s ) + 2a 1 ln(s/s ), 



(2) 



with a' « 0.25 GeV~ 2 . Thus, the radius of strong inter- 
actions is expected to be a factor of 1.5 larger at LHC as 
compared to fixed target energies. 

In the partonic picture, the mechanism for the increase 
of the radius of strong interactions with energy is the so- 
called Gribov diffusion. The emission processes in the 
soft parton ladder give rise to a random walk of partons 
in the transverse plane, reminiscent of a diffusion pro- 
cess [8| . If one writes the amplitude of elastic scattering 
of two hadrons as a product of two i-dependent form fac- 
tors, each parameterized by the transverse radius of the 
hadron, R, 
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exp(i.R 2 /4) exp(ti^/4), 



(3) 



one can interpret the shrinkage of the diffractive cone as 
being due to an increase of the transverse spread of par- 
tons. In terms of the average parton momentum fraction, 
x, this implies that the transverse area occupied by the 
low virtuality partons in the hadron, R 2 , increases with 
decreasing x roughly as 



R 2 (x) = Rl + 2a'ln(x /x) 



(4) 



Detailed information about the distribution of the pp 
cross section (both elastic and inelastic) over impact pa- 
rameters can be obtained from the impact parameter rep- 
resentation of the pp elastic scattering amplitude, see e.g. 



Ref. 0]. We write the invariant elastic amplitude in the 
form 

A pp (s,t) = — [ d 2 be~^ ±b) T pp (s,b) (5) 

47T J 

= ~2 I dbbJ o(^±b)TPP(s,b), (6) 

where A_l is a transverse momentum vector, with t — 
— Aj_ and Aj_ = |Aj_|, and Jo denotes the Bessel func- 
tion. Our normalization of the amplitude is the same 
as in Ref. 0, c/. the relation to the total cross section, 
Eq. below. The dimensionless complex function T pp 
is called the profile function of the elastic amplitude. In 
this representation the integrated cross section for elastic 
scattering is given by 



dt \A pp (s,t)\' 



d 2 b \T pp (s,b)\ 2 



(7) 
(8) 



In the last step, we have used that at large s the lower 
limit of the i-integral can be replaced by — oo, and that 
for an amplitude independent of the azimuthal scattering 
angle 



dt = 4tt 



d 2 A A 



(9) 



A similar representation can be derived for the total cross 
section for pp scattering, which by the optical theorem is 
proportional to the imaginary part of the forward (t = 0) 
elastic amplitude: 

a pp t {s) = -ImA pp (s,t = 0) (10) 
s 

= 2 J d 2 bReT pp (s,b). (11) 

Finally, taking the difference of Eqs. i|ll[l and (JSJ), one ob- 
tains a representation of the inelastic (total minus elastic) 
pp cross section as an integral over impact parameters: 

e, af£( S )-<( S ) 

= J d 2 b [2ReT pp {s,b)-\r pp (s,b)\ 2 ] . (12) 

The integrand of Eq. (|12fl represents the distribution 
of the cross section for generic inelastic collisions (i.e., 
summed over all inelastic final states) over impact pa- 
rameters. It is convenient to define a normalized b— 
distribution as 



Pin(s,b) 



2Re T pp (s,b) - \TPP{s,b)\ 2 



(13) 



For a quantitative estimate of this distribution we can 
use phenomenological parameterizations of the pp elastic 
scattering amplitude, which fit the presently available pp 
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FIG. 3: The normalized impact parameter distribution for 
generic inelastic collisions, Pi n (s,b), Eq. 113L obtained with 
the parameterization of the elastic pp amplitude of Islam et 
al. |t| ( "diffractive" part only). The plot shows the "radial" 
distribution in the impact parameter plane, 27r6Pi n (s, 6). The 
energies are = 500 GeV (RHIC) and 14000 GeV (LHC). 



elastic data at collider energies. The results obtained 
with the parameterization of of Islam et al. ( "diffrac- 
tive" part only) are presented in Fig. 01 for energies cor- 
responding to RHIC and the LHC. It should be noted 
that the predictions for LHC are based on extrapolation 
of fits to the presently available data over nearly two or- 
ders of magnitude in s. The biggest uncertainty in the 
extrapolation appears to be due to the uncertainties in 
the measurement of a^ t at the Tevatron and the limited 
range of t covered in the collider measurements of elastic 
pp scattering. 

In principle one should include here also effects of in- 
elastic diffraction. However, this contribution should be 
rather small at LHC energies due to the blackening of in- 
teraction, see Ref. ^3 f° r a recent review. Besides this, 
one expects that a significant part of inelastic diffraction 
at t < is due to the spin flip amplitudes (llj| . 



III. TRANSVERSE SPATIAL DISTRIBUTION 
OF HARD PARTONS IN THE NUCLEON 

We shall now review what is known about the trans- 
verse distribution of hard partons in the nucleon. Based 
on this, we shall later proceed to estimate the impact 
parameter distribution in pp collisions with hard dijet 
production. 

Numerous measurements of hard inclusive scattering 
processes (DIS, Drell-Yan pair production) have pro- 
duced a rather detailed picture of the longitudinal mo- 



mentum distribution of partons in the nucleon. The 
study of the transverse spatial distribution of partons is 
still at a much more primitive stage. Information about 
the transverse spatial distribution of gluons is contained 
in the so-called two-gluon form factor of the nucleon, 
which parametrizes the ^-dependence of the (generalized) 
gluon distribution in the nucleon, 



?0M;Q 2 ) = g(x;Q 2 ) F g (x,t;Q 2 ), 



where 



F g (x,t = 0;Q 2 ) 



1. 



(14) 



(15) 



and g(x; Q 2 ) is the usual gluon distribution in the nu- 
cleon. We define the Fourier transform of this form factor 



F g (x,p; Q 2 ) 



d 2 A ± 



'< A ^> F g (x,t = -Ai;Q 2 ), 

(16) 

where p is a two-dimensional coordinate variable, and 
p = \p\ (this variable is named b in Ref. [l^; m this paper 
we reserve b for the impact parameter vector of the pp 
collision) . This function describes the spatial distribution 
of gluons in the transverse plane. It is normalized to unit 
integral over the transverse plane, 



J d 2 pF g (x, P -Q 2 ) = 1. 



(17) 



A measure of the transverse size of the nucleon for given 
x and Q 2 is the average of p 2 calculated with this dis- 
tribution, which is identical to 4 times the i-slope of the 
two-gluon form factor at t = 0, 



(p 2 }(x,Q 2 ) = / d 2 pp 2 F g {x, P] Q 2 ) 







4 -F g (x lt ;Q 2 ) 



(18) 
(19) 



For sufficiently small x (< 0.3) the parameter p can be 
interpreted as the distance of the parton from the center 
of mass of the nucleon in the transverse plane. For larger 
x the interpretation of the p-distribution becomes less 
intuitive, as in this case one can no longer neglect the dif- 
ference between the longitudinal momentum of the spec- 
tator system and that of the whole nucleon. In the limit 
x — ► 1 the active parton would carry the entire longitu- 
dinal momentum of the nucleon, and only soft partons 
would be left in the spectator system, see Ref. [L3] for a 
discussion. This shall not concern us here, since we shall 
be interested in the gluon distribution at x < 0.05, as 
relevant for hard dijet production in the central rapidity 
region (see below). 

At moderately small x (> 0.001), it is possible to ob- 
tain information about the two-gluon form factor at a 
resolution scale of Qq ~ 2 4 GeV 2 from the analysis of 
exclusive J/tjj photo (or electro-) production off the nu- 
cleon It turns out that for x > 0.01 this form factor is 
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FIG. 4: Our model for the ^-dependence of the average trans- 
verse gluonic size squared of the nucleon, (p 2 ) at the scale 
Qo — 2 -T- 4GeV 2 relevant to J/t/j production. Short-dashed 
line: {p 2 } = 0.28 fm 2 , as extracted from the f-slope of the Jhb 



production cross section measured in various experiments _ 
Long-dashed line: Sum of the constant value (p 2 ) = 0.28 fm 2 
and the pion cloud contribution calculated in Ref. [l2Tl . Solid 
line: The parameterization Eq. 12611 . based on the experimen- 
tal value of ah ard , Eq. 125H 15(. 



significantly harder than the electromagnetic form factor 
of the nucleon. The i-dependence of the cross section is 
well described by a dipole form factor 



F g (x,t) 



with a mass parameter 



{1-t/mir 



(20) 



1.1 GeV 2 > mt 



(a; > 0.01). (21) 



The corresponding spatial distribution of gluons in the 
transverse plane, Eq. I|16|) . is given by 

F g (x,p) = ^ (!?f) Kl (m g p), (22) 

where Ki denotes the modified Bessel function. Note 
that this function is positive, in agreement with the gen- 
eral positivity condition for the transverse coordinate- 
dependent gluon density, derived in Ref.jl3|. The aver- 
age (p 2 ) corresponding to this distribution is inversely 
proportional to the mass parameter squared, 



(P 2 ) = 4" 
g 

Numerically, the value (|21|l amounts to 

(p 2 ) w 0.28 fm 2 (x > 0.01), 



(23) 



(24) 



which is a factor of ~ 1.5 smaller than the proton's 
"transverse electric charge radius squared", (2/3)(r 2 ) om . 

For smaller values of x (0.01 < x < 0.1) the transverse 
gluonic size of the nucleon starts to increase. This can 
be explained semi-quantitatively by the kicking in of con- 
tributions from the pion cloud of the nucleon, which are 
suppressed for x > M^/Mm [l2|, see Fig. The analy- 
sis of J/ip photoproduction data has shown that the size 
keeps growing also for x < 0.01 ^ n au > tne rate °f 

the increase of the gluonic size between x ~ 10~ 2 and 
10~ 3 is about a factor of two smaller than that of the 
total cross section, which is dominated by soft physics 
(see Section [nj: 



1 d{ P 2 ) 
4 d\n{l/x) 



a hard ? 
(x = 10 



* 0.125 GeV" 2 

3 -io- 2 ), 



(25) 



which should be compared with Eq. J2J) . We can parame- 
trize the ^-dependence of (p 2 ), at the scale Q\ probed in 
J/?/ 1 production, by combining this experimentally de- 
termined rate of increase with our model estimate of 
(p 2 ) (constant plus pion cloud contribution) at x = 0.1, 
(p 2 ) = 0.31 fm 2 12]. This amounts to the parametriza- 
tion 



0.31 fm 2 
0.28 fm 2 



0.0194 fm 2 ln%5- 



(26) 

This simple form fits well the x-dependence of (p 2 ) due to 
pion cloud contributions in the region 0.01 < x < 0.1, as 
calculated in Ref. [13, see Fig.il and continues it down 
to smaller values of x using the experimentally measured 
rate of increase. 

For our estimates of the probability of hard multijet 
production we need to model not only the average trans- 
verse size of the nucleon, but the full transverse spatial 
distribution of gluons. For simplicity, we shall assume 
that at the scale Q\ the p-distribution at all relevant x 
can be described by the Fourier transform of a dipole 
form factor, cf. Eq. I|22|) . but with an x-dependent mass 
parameter. This parameter is then uniquely determined 
by the value of (p 2 ), as given by Eq. via Eq. 

This defines our model of the x- and p-dependent distri- 
bution of gluons at the scale Qq. 

We are interested in the p-dependent gluon distribu- 
tion at large virtualities, corresponding to hard dijet pro- 
duction at LHC. This requires to take into account the 
effect of DGLAP evolution on the p-dependent distribu- 
tions. The Q 2 -evolution of the parton distributions is 
diagonal in p. It degrades the longitudinal momentum 
fractions x of the partons, resolving a parton with given 
x into a collection of partons with smaller x, while the 
transverse location of the new partons is practically the 
same as that of the "parent" parton provided Q ^> 1/p. 
Nevertheless, Q 2 evolution does change the p-profile of 
the distributions at given x. Generally speaking, evolu- 
tion will reduce the rate of broadening of the distributions 
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FIG. 5: The change of the normalized p-profile of the FIG. 6: The change of the average transverse gluonic size 



gluon distribution, F g (x, p;Q 2 ), Eq. 1161 . with Q 2 , as due 
to DGLAP evolution, for x — 1CF 3 . The input gluon dis- 
tribution is the Gluck-Reya-Vogt parametrization at Qq = 
3 GeV 2 , with a dipole-type p-profile, Eq. 12211 . of size deter- 
mined by the parametrization Eq. i'Aiil . 



with decreasing of x. This happens because with increas- 
ing Q 2 the parton distributions at the high scale become 
sensitive to the input distribution at the low scale Qq 
at larger and larger x-values, where their transverse size 
becomes small. 

We have studied numerically the leading-order Q 2 - 
evolution of our model of the p-dependent gluon distribu- 
tion at the scale Qq cf. Eqs. (J23), i(2T)|) and employ- 
ing the numerical method described in Ref . • We have 
used the Gliick-Reya-Vogt leading-order parameteriza- 
tion [17| to describe the total (integrated over p) gluon 
and singlet quark distributions at the input scale Qq and 
modeled their p-profile as described above. The value of 
the input scale we have taken as Qq = 3 GeV 2 , which is 
the central value of the range of scales associated with 
J/ip photoproduction. For simplicity we assume identi- 
cal p-profiles for the gluon and singlet quark distribu- 
tions at the input scale. The results of the DGLAP evo- 
lution of the p-dependent gluon distribution are shown 
in Figs. and Fig. Fig. [3] shows the normalized p- 
profile, F g (x, p;Q 2 ), cf. Eq. (|16fl . corresponding to the 
distribution after evolving to the higher scale Q 2 , for the 
value x = 10~ 3 . Fig. [|)]shows the Q 2 -dependence of the 
average size squared, (p 2 ), as induced by DGLAP evolu- 
tion, for values of x = 10~ 2 , 10~ 3 and 10 -4 . The corre- 
sponding curves for the singlet quark distribution would 
be qualitatively similar. One can see that the effect of 
transverse broadening of the gluon distribution with de- 
creasing x, which is rather small already at the initial 
resolution scale, is further reduced at the scale relevant 



squared, (p }, due to DGLAP evolution, for x — 10 
and 1CT 4 . 



10' 



for the LHC kinematics {Q > 20 GeV). 

Fig.|Sldemonstrates that the deviations of the p-profile 
from the dipole shape at the initial scale, Eqs. (|22|l . due 
to DGLAP evolution are very small. This suggests a 
simplified parametrization of the combined Q 2 - and x- 
dependence of the p-profile of the gluon distribution, in 
which the dipole shape is assumed to hold at all Q 2 (and 
x), and the Q 2 - (and x-) dependence is entirely ascribed 
to the mass parameter m 2 . We fit the combined Q 2 - 
and ^-dependence of (p 2 ) due to DGLAP evolution in 
the region Qq<Q 2 < 10 6 GeV 2 , as shown in Fig. HO by 
the simple two-parameter form 



(27) 



where Q 2 , = 3 GeV 2 , (p 2 )(x, Q 2 ,) is defined by Eq. ®, 
and 



A = 1.5, 



1 



a = 0.0090 In-. (28) 

x 



For each x and Q 2 , this value of (p 2 ) defines a dipole 
mass parameter m 2 via Eq. 1)23(1. Our model for the p- 
dependent gluon distribution is then given by the Gliick- 
Reya-Vogt leading-order parametrization for the total 
distribution at the scale Q 2 , times the normalized dipole 
p-profile, Eq. (|22|l . with this mass parameter. This 
parametrization has the correct Q 2 -dependence of the 
transverse size is "built in" , removing the need to perform 
explicit DGLAP evolution of the p-dependent distribu- 
tions, which is very convenient for the following studies. 
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IV. IMPACT PARAMETER DISTRIBUTION 
FOR A HARD MULTIJET TRIGGER 

Using the information about the transverse distribu- 
tion of partons in the proton, we can now investigate the 
impact parameter dependence of the cross section for in- 
elastic collisions with production of two jets in the central 
rapidity region. In particular, we shall show that a trig- 
ger on hard dijet events allows to reduce the effective 
impact parameter of pp collisions as compared to generic 
inelastic collisions studied in Section ITT1 

We consider the production of two jets (with equal 
but opposite transverse momentum) in a binary parton- 
parton collision. The resolution scale is given by the 
transverse momentum squared of one the jets, q\. The 
momentum fractions of the two colliding partons with 
respect to their parent protons, x\ and X2, can be recon- 
structed from the measured energy and momenta of the 
two jets. Four-momentum conservation implies for the 
scattering at 90° at the center of mass of two partons 



Xi x 2 



4qj 

S 



(29) 



where Aq^_ is the invariant mass squared of the two-jet 
system. In the following we shall be interested in two jets 
near zero total rapidity, which requires 



■r-2 ■ 



(30) 



Under this condition the momentum fractions are com- 
pletely fixed by Eq. I|29|l . In the following we consider the 
dijet production due to collision of two gluons, since such 
partonic collisions give the dominant contribution to the 
total cross section. The probability for a binary collision 
of two gluons is proportional to the product of the gluon 
densities in transverse space in the two colliding protons, 
taking into account that their transverse centers are sep- 
arated by a distance b — the impact parameter of the 
pp collision, see Fig. [7] This implies that the distribu- 
tion of the cross section for such events over the impact 
parameter b is given by 



ft (6) 



Pi 



d 2 P2 8&(b- 



Pi + Pi) 



X Fg{Xl,px) Fg(Xl,p 2 ), 



(31) 



where x\ = 2q±/^/s, cf. Eqs. (|2*9"jl and J3UJ, and the scale 
of the gluon p-profiles is q\ . This distribution is normal- 
ized such that the integral over all b is unity. Since it has 
the form of a convolution in the parton transverse posi- 
tions, it can also be expressed as the Fourier transform 
of the square of the two-gluon form factor, Eq. (|20|) . In 
particular, for the two-gluon form factor of dipole form, 
Eq. <E2), used in our model of the p-dependent gluon 
distribution (see Section HTT|l one obtains 



Lztt 



K 3 (m g b), 



(32) 




FIG. 7: Illustration of the overlap integral of parton distribu- 
tions in the transverse plane, defining the ^-distribution for 
binary parton collisions producing a dijet, Eq. 131> . 
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FIG. 8: The b-distribution for the trigger on hard dijet pro- 
duction, Pi{b), obtained with the dipole form of the gluon 
^-profile, Eq. for ^ = 14000 GeV and q± = 10 GeV 

and 100 GeV. The plots show the "radial" distributions in 
the impact parameter plane, 2nbP2(b). Also shown is the 
corresponding distribution for a trigger on double dijet pro- 
duction, Pa{V), with the same p±. 



where m g should be substituted by the value correspond- 
ing to x\ = 2qj_/^/s and Q 2 = q\, see Eq. l|2*7jl . 

Fig.[S]shows the distribution P2{b) for a center-of-mass 
energy of = 14000 GeV (LHC), and two values of the 
jet momentum, q± = 10 GeV and 100 GeV. One sees 
that the distribution is rather insensitive to the precise 
value of the jet momentum. This can be explained by 
the relatively slow decrease of (p 2 ) with increasing x and 
Q 2 . The average values of impact parameter squared, 
(b 2 ), calculated with these distributions, is 0.71 fm 2 for 
q±_ = 10 GeV and 0.63 fm 2 for q ± = 100 GeV. 

In Fig. |31 we assume production of a two-jet system 
at zero rapidity, cf. Eq. I|30|l . If we considered instead a 
two-jet system at some non-zero rapidity, y, the (anyway 
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weak) dependence of the ^-distributions in Eq. i|31|l on x\ 
and X2 would work in opposite directions, leading to an 
extremely weak dependence of our results on the rapidity 
of the produced system over a wide range of y. 

In Fig. we compare the 6-distribution for the hard 
dijet trigger, P2{b) (solid line), with the 6-distribution for 
generic inelastic events, P; n (s), estimated in Section ITU 
The short-dashed line in Fig. [5| shows the distribution 
Pi n (s) obtained from the parameterization of the elas- 
tic pp amplitude of Islam et al. 9} ( "diffractive" part 
only). Shown are the results for y^i = 14000 GeV (LHC), 
1800 GeV (Tevatron pp), and 500 GeV (RHIC). A mo- 
mentum of q± = 25 GeV was assumed for the dijet trig- 
ger. One sees that in all cases the 6-distribution for dijet 
events is much narrower than the one for generic inelas- 
tic collisions. The corresponding averages (b 2 ) are given 
in Table [TJ The average (b 2 ) for the hard dijet trigger 
rises much more slowly with s than for generic inelastic 
collisions, which are dominated by soft physics. Thus, 
the reduction in effective impact parameters due to the 
dijet trigger is most pronounced at LHC energies, where 
(b 2 ) is reduced to ~ 1/4 its value for generic inelastic 
collisions. 

A further reduction of the effective impact parameters 
can be achieved by a trigger on events with two dijets, 
i.e., two binary hard parton collisions (such processes 
can be easily distinguished from the leading twist 2^4 
processes in the collider experiments, see e.g. Ref. 0). 
It was estimated in Ref. [l£| that this reduces (fe 2 ) by 
a factor of two as compared to the single dijet trigger. 
In our approach, the ^-distribution for the double dijet 
trigger is given by 



P 4 (6) 



Pi(b) 



d 2 bp 2 (b) 



(33) 



For simplicity we assume here identical x\ and q± for the 
two dijets; the definition could easily be generalized to 
allow for different values. For the two-gluon form factor 
of dipolc form, Eq. flffi. this becomes 



P 4 (6) 



7 m 2 
"367 



rrigb 



[K 3 (m g b)Y 



(34) 



where again m g should be substituted by the value cor- 
responding to X\ = 2q±/^s and Q 2 = q\, see Eq. 
Fig. |S] shows that the ^-distribution for the double dijet 
trigger, P±(b) is equally insensitive to the precise value of 
q± as that for the dijet trigger, P2(b). The comparison 
in Fig. [U] and Table [I] shows that the double dijet trig- 
ger allows for a further reduction of the effective impact 
parameters by a factor of ~ 2.5 compared to the dijet 
trigger. 

In calculating Pi(b) we have made the assumption that 
the gluons are not correlated in the transverse plane. To 
test this assumption we can compare the rate of double 
dijet production in our model to the one which is ob- 
served in the CDF experiment • The ratio of the cross 



Facility 


^5/GeV 


(& 2 ) 2 /fm 2 


<& 2 }4/fm 2 


(& 2 ) in /fm 2 


LHC 


14000 


0.67 


0.26 


2.7 


Tevatron 


1800 


0.63 


0.24 


1.8 


RHIC 


500 


0.59 


0.23 


1.43 



TABLE I: The average impact parameter squared, (b 2 ), cor- 
responding to the fa-distributions P-2(b),P4(b), and P m (&), 
shown in Fig. [§] 



section of double dijet events and the square of the single 
dijet cross section is proportional to 0^1 



Ccff 



d 2 bP 2 (b) 



(35) 



In our calculation we find a D g = 34 mb, which should be 
compared with a c s = 14.5 ± l.Ttgijjnib reported by CDF 
assuming that there is no correlations in the longitudinal 
distribution of partons. (In Ref. [19| cr c ff ~ 30 mb was ob- 
tained assuming that the parton distribution similar to 
that of valence quarks, a hypothesis resembling the con- 
clusion we derived from the leading twist analysis of the 
J/ip elastic photoproduction 3].) A factor of two differ- 
ence between the theoretical number and the data may 
indicate that there are significant transverse correlations 
in the parton density at the resolution scale of Q > 5 GeV 
probed by CDF. Such correlations could result, for exam- 
ple, due to the DGLAP evolution from a low Q 2 scale of a 
couple GeV 2 to Q ~ 25 GeV 2 , since the partons emitted 
in the course of this evolution would have small trans- 
verse separation (the so-called "hot spots" of Ref. [20j|L 
Assuming that the difference between the uncorrelated 
model and the data is due to such local correlations in 
b we would obtain a 6-distribution for the double dijet 
trigger approximately as 



Pi.corr(b) W P 2 (&) 

+ Pi{b) 



cr Gff (model) - cr cff (CDF) 

(T ff (model) 
CTcff(CDF) 



(T c ff (model) 



(36) 



where P2Q)) and Pi(&) are the above uncorrelated model 
estimates. It is clear from the inspection of Fig. that 
this still corresponds to a large reduction (by a factor 
~ 1.5) of the effective impact parameters with the double 
dijet trigger as compared to the single dijet trigger. 

It should be noted that the jet cross sections discussed 
here are sensitive to the choice of "primordial" trans- 
verse momentum distribution of the colliding partons, to 
higher-order radiative corrections, and to the jet defi- 
nition adopted. To our knowledge, the same procedure 
was used in analysing single and double dijet events in 
Ref. (one of the jets was taken to be a photon for 
simplicity), so most of these effects should cancel in the 
cross section ratio. An additional source of uncertainty is 
the correlation of the soft background in hadron produc- 
tion with the centrality of the pp event, as discussed in 
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FIG. 9: Solid lines: fa-distributions for the dijet trigger, -p2(&), with q± = 25GeV, as obtained from the dipole-type gluon 
p-profile, Eq. 13211 . Long-dashed line: 6-distribution for double dijet events, PiQj). Short-dashed line: ^-distribution for generic 
inelastic collisions, obtained from the parameterization of the elastic pp amplitude of Islam et al. ( "diffractive" part only), 
cf. Fig. 01 Shown are the results for Vs = 14000 GeV (LHC), 1800 GeV (Tevatron pp), and 500 GeV (RHIC). The plots show 
the "radial" distributions in the impact parameter plane, 2nbP(b). 



Section I VII A quantitative study of parton-parton cor- 
relations certainly requires a more careful investigation 
of these effects. 

To summarize, we have demonstrated that a trigger 
on events with (one or more) dijets near zero rapid- 
ity strongly reduces the effective impact parameters in 
pp collisions at LHC energies. Such a trigger can thus 
be used as a "centrality filter" . This is of considerable 
practical interest, as the characteristics of the final state 
strongly depend on the centrality of the pp collision. 



V. APPROACHING THE BLACK-BODY LIMIT 
IN CENTRAL PP COLLISIONS 



An interesting feature of central pp collisions at high 
energies is that large-x partons (x > 0.01) in one nucleon 
pass through a strong gluon field in the other nucleon. 
This field can become so strong that the interaction of the 
parton with the other nucleon approaches the Black Body 
Limit (BBL), in which the probability for inelastic scat- 
tering becomes unity, and the cross section becomes com- 
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parable to the transverse size of the strong gluon field. 
This phenomenon would have dramatic consequences for 
particle production in the forward region (dilepton pro- 
duction, hadron multiplicities and transverse momentum 
distributions), which have become the subject of inten- 
sive theoretical investigation. In this section we want to 
quantify the proximity to the BBL for pp collisions at 
LHC energies. Specifically, we want to show how a trig- 
ger on hard dijet (or multi-jet) production, which reduces 
the effective impact parameters in pp collisions, greatly 
increases the region (in the momentum fraction, x, and 
the virtuality, Q 2 ) in which partons experience interac- 
tions close to the BBL. For a recent review of approaches 
to "taming" the growth of parton densities at small x, 
based on the impact parameter eikonal approximation 
and the leading-loga; approximation, see Refs. pl| . 

To simplify the discussion, we consider instead of the 
scattering of a colored parton the scattering of a small 
color-singlet dipole off the other nucleon. This is in the 
spirit of the dipole picture of high-energy scattering of 
Mueller [2^J . The distribution of the inelastic and elas- 
tic cross sections, af p and <r^, over the dipole-proton 
impact parameter, p, can be expressed in terms of the 
profile function of the dipole-nucleon elastic scattering 
amplitude, T d P(s,p) (cf. Section HU: 



dp 



in. cl 



00 



d2 P a in.el( S iP)i 



with 



v$(s,p) = 2ReT d P( S ,p)-\T dp ( S ,p)\ 2 



(37) 



(38) 
(39) 



The functions a-J^ , (s, p) are dimensionless and can be 
interpreted as the "cross sections per unit transverse 
area". If the target proton were a "black" disk of ra- 
dius R, the probability for a dipole hitting the target 
to undergo inelastic scattering would be unity, and one 
would have of p (s,p) = 1 for p < R. This would imply 



T dp {s,p) = 1 for p<R. 



(40) 



From Eq. I|38l) it follows that in this case the elastic cross 
section would also be unity 



a ei ( s > P) = 1 for p < R, 



(41) 



and thus the total (elastic plus inelastic) cross section 
would be twice the inelastic one. For realistic inter- 
actions, one may say that the proton becomes black if 
T dp (s,p) approaches unity in a certain region of impact 
parameters. In the following we shall use this criterion 
to quantify the approach of dipole-nucleon interactions 
to the BBL. A similar approach was used in Ref. |2^| to 
study the proximity of 7* AT scattering to the BBL. 

For small dipole sizes the cross section for inelastic 
dipole-nucleon scattering at fixed impact parameter b is 



riven by the leading-twist perturbative QCD expression 
11 



Ctt 2 



LTV" j r i 3 

Here C is a color factor, 

C = 



d 2 a s (Q 2 )x'g(x',p;Q 2 ). (42) 



1, 33 dipole 
9/4, 88 dipole 



(43) 



d is the dipole size, a s (Q 2 ) the leading-order QCD run- 
ning coupling, 



a s (Q 2 ) 



47T 



Po \n(Q 2 /A 2 QCD Y 
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~-N, 



-N 



(44) 
(45) 



(Aqcd is the scale parameter), and g(x, p; Q 2 ) the impact 
parameter-dependent leading twist gluon density in the 
nucleon, 



g{x' lP -Q 2 ) = 5 (x',Q 2 ) F g (x',p;Q 2 ), 



(46) 



where g(x',Q 2 ) is the usual (total) gluon density and 
F g (x' , p;Q 2 ) the normalized p-profile, Eq. (|Pd|| . The 
scale parameter Q 2 in a s and the gluon density is re- 
lated to the dipole size by 



I 2 ' 



(47) 



where A is a dimensionless parameter whose value is to 
be determined from phenomenological considerations. A 
value of A = 9 sets d equal to the average dipole size con- 
tributing to the longitudinal photon-nucleon cross sec- 
tion at large Q 2 ; we shall use this value of A in the follow- 
ing. The gluon momentum fraction probed in inelastic 
dipole-nucleon scattering is determined by the invariant 
mass squared of the produced system, M 2 , 



M 2 + Q 2 



(48) 



Here we have in mind generic inelastic production with 
M 2 ~ Q 2 , so we take 



2Q 2 



= 2x 



(49) 



as the argument for the gluon distribution. In the fol- 
lowing we shall regard the cross section viz. the profile 
function as a function of x = Q 2 / s and the dipole-proton 
impact parameter, p. 

Comparing Eqs. 139|) and (|42|l we obtain an equation 
for the profile function T (regarded now as a function of 
x and p) corresponding to the leading-twist QCD result. 
If we assume the elastic amplitude to be imaginary, i.e., 
r to be real, we obtain 



2 Y dp — (r d P) 2 = (j d P 



in, lt- 



(50) 
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FIG. 10: The profile function for elastic dipole-nucleon scattering, V dp (x,p), obtained from Eq. I|51fl with the inelastic cross 
section given by the leading-twist expression, Eq. 1421 . Shown are the results for an 88 dipole with Q 2 = 100 GeV 2 (upper left 
panel), 50 GeV 2 (upper right panel), and 20 GeV 2 (lower left panel), for various values of x. 



The relevant solution of this quadratic equation is 



T*>(x,p) = l-yJl-cr£ LT (x,p). (51) 

Note that this solution implies that T dp — > 0, and thus 
Coi — > [cf. Eq. ffiKjl]. for a- m , lt — > 0, i.e., in the limit 
of small dipole size. For the corresponding solution with 
the "+" sign in front of the square root the elastic cross 
section would tend to a constant value for zero dipole 
size. 

We evaluate the leading-twist expression for the inelas- 
tic 88 dipole-nucleon cross section (|42|) with the model 



for the p-dependent gluon distribution described in Sec- 
tion IIIII based on the dipole form factor with x- and 
Q 2 -dependent mass parameter. Fig. ^| shows the pro- 
file function T dp (x,p) obtained from Eq. I|51|) for values 
x = 1CT 2 , 1(T 3 and 1(T 4 , and Q 2 = 100, 50 and 20 GeV 2 ; 
the dipole size in each case is determined as ar — X/Q 2 
with A = 9. The results clearly show the approach to the 
BBL for decreasing x (as a result of the growth of the 
gluon density) and/or decreasing Q 2 (as a result of the in- 
creasing dipole size) . Note that the results for the profile 
function shown here are physically meaningful only in the 
region of impact parameters where T dp is significantly less 
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FIG. 11: The critical value of impact parameter, p C rit, for 
which the profile function of elastic dipole-nucleon scattering, 
V dp {x, p), exceeds the value F^ it = 0.5, cf. Eq. Here p C rit 

is shown as a function of x. A value of p crr t = implies that 
T dp (x,p) < T d f it for all values of p. Shown are the results for 
an 88 dipole; the Q 2 -values are the same as in Fig. I1UI 



than unity; for T dp ~ 1 the simple leading-twist formula 
for the inelastic dipole-nucleon cross section, Eq. (|42[1 . 
breaks down. In other words, the figures should be read 
as indicating the region in impact parameter for which 
the cross section approaches the BBL, rather than how 
it is approached once the interaction is sizable. An ob- 
jective measure of the size of this region is the value of p 
at which F dp exceeds a certain critical value, r crrt , 



r dp (x, P ) > r 



dp 



for 



P < Pcrit- 



(52) 



Fig. 1111 shows p cr it as a function of x, and for the above 
values of Q 2 , for r* it = 0.5. This value of r dp corre- 
sponds to a probability of having no inelastic interaction 
of only 0.25, cf. Eq. I|39f) . which is a reasonable guideline. 
In Fig. ^]a value of p cr i t = implies that T dp (x,p) is 




FIG. 12: The probability for partons in one proton with given 
x and Q 2 to interact with the other proton near the BBL, as 
defined by Eq. (I53II . as a function of the impact parameter 
of the proton-proton collision, b. Shown are the results for 



smaller than T c f it for all values of p; even in the center of x = 10 (upper panel) and x = 10 (lower panel), for 
the nucleon where the gluon density is maximum (p = 0). Q — 10, 20 and 30 GeV . 



We now turn to the question under which conditions 
in proton-proton collisions a parton in one proton will 
interact with the gluon field in the other proton near 
the BBL. The probability for this to happen depends 
on the transverse position of the parton relative to the 
center of its parent proton, pi, as well as the impact 
parameter of the proton-proton collision. An interesting 
measure is the total probability for partons in one proton 
with given momentum fraction x and virtuality Q 2 (but 
arbitrary transverse position) to interact with the other 
proton near the BBL. It is given by the overlap integral 



of the normalized transverse spatial distribution of the 
partons in their parent proton (shifted by the impact 
parameter vector, 6), with the characteristic function of 
the "black" region in the other proton, 0(pi < Pcrit), 

fbiack(6) = J d 2 pi 8(pi < pent) F g (x,p 2 ) 

(pa = \px-b\). (53) 

This integral measures the fraction of the partons with 
given x and Q 2 which hit the other proton in the "black" 
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central region. Here p cr i t depends on x and Q 2 , see above. 
Generally, the probability Pbiack(^) is maximum for cen- 
tral collisions (b — 0) and decreases with increasing b. 
Note that if the target proton were a "black" disk of ra- 
dius R, and the parent proton of the parton a disk of 
same size, Pbiack would be unity at b — 0. Fig. 1121 shows 
-Pbiack(fr) for partons with virtuality Q 2 = 20GeV 2 and 
different values of x. 

In actual hadron-hadron collisions, the parton which 
probes the gluon field in the other proton is resolved by 
a hard collision with a parton in the other proton, re- 
sulting in hadron production, see Fig. |2 Let xr be the 
momentum fraction of the resolving parton. Production 
of hadrons with transverse momentum p± then resolves 
partons with momentum fraction [cf. Eq. (f3~|) ] 



and virtuality 



XRS 



q 2 = m 



(54) 



(55) 



We are interested in partons with relatively large xr (~ 
10 _1 ), which are able to resolve small-cc partons in the 
other proton whose interactions are close to the BBL. 
Generally speaking, large p± select partons with large x 
and large Q 2 , for which interactions close to the BBL 
are unlikely Thus, the probability of interactions close 
to the BBL decreases with increasing p±. On the other 
hand, a certain minimum value of the parton virtuality, 
Q 2 , and thus of p±, is required for the concept of parton 
resolution by jet production to be applicable. Thus, it 
is crucial to establish that there is a "window" in p± in 
which one is sensitive to BBL effects while at the same 
time the partonic description is still applicable. 

An important quantity is the maximum value of p± 
(for given xr) for which the resolved parton sees the 
other proton as "black" . We can estimate this maximum 
p± with the help of the probability Pbiack introduced in 
Eq. H53JI . For given xr, and given impact parameter b, 
we ask for the maximum value of p\ for which -Pbiack(fr) 
exceeds a certain critical value: 



fblackO) > Pm 



for p\ < pIbbl- 



(56) 



Fi gure 1131 shows p 2 ^ BB l corresponding to the criterion 
that -Pbiack(&) > 0.5, as a of the impact parameter of the 
proton-proton collision, b. One sees that the maximum 
p± drops rapidly with increasing b. [For b > 1 fm the 
maximum value of p\ defined according to Eq. (|56|l be- 
comes < 1 GeV 2 and thus physically meaningless.] For 
central collision, b < 0.5 fm, the values of p 2 ^ bbl are 
larger than 10 GeV 2 , so that the underlying leading-twist 
approximation is well justified. For such impact param- 
eters it is possible to explore the properties of the BBL 
using resolved hard partons as a well-defined probe. 

Strictly speaking, the scattering of the small-x parton 
in one proton with the large-a^ parton in the other pro- 
ton occurs at finite cm. angles, and thus leads to a loss 




FIG. 13: The maximum value of the jet transverse momentum 
squared, p\ BBL , for which a small-x parton in one proton, 
resolved in a collision with a large-a;^ parton in the other 
proton, interacts with the other proton close to the black- 
body limit (BBL), as a function of the impact parameter of 
the proton-proton collision, b. The criterion for proximity 
to the BBL is Pbiack(fr) > 1/2, see Eq. 15611 . Shown are the 
results for xr = 0.1, 0.2 and 0.3. 



of light-cone fraction for the large-rr^ parton. This effect 
should be included in a more accurate treatment. 

The above results clearly show that in order to explore 
the BBL with partons resolved in backward/forward jets 
as a probe one needs to limit the effective impact parame- 
ters in the pp collisions. Here the idea of using hard dijet 
production at central rapidities as a "centrality filter" , 
developed in Sect ion llVl comes in. Assuming that the si- 
multaneous production of central and forward/backward 
jets can be described incoherently, i.e., in a probabilistic 
manner, the probability for interactions near the BBL in 
events with hard dijet production is given by the average 
of the impact parameter-dependent p_l,bbl distribution 
of Fig. E| with the effective 6-distribution implied by the 
hard dijet trigger, P2(b), see Section ITVl 



(Pj_,bbl>2 



J d 2 bp\ 



BBL 



(b) P 2 (b). (57) 



A similar definition applies to the double dijet trig- 
ger with 6-distribution Pi(b). For the LHC energy, 
y/s = 14000 GeV, and a dijet trigger with momentum 
p± = 25 GeV, cf. Fig. El the average values of p\ BBL 
obtained from Fig. ^| are shown in Fig. ^] (solid line) as 
a function of the resolving parton's momentum fraction, 
xr. Also shown are the corresponding averages with a 
double dijet trigger with the same p± (dashed line). One 
sees that the average values of p\ BBL are all ^> 1 GeV 2 , 
i.e., in the region where our assumption of resolved hard 
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FIG. 14: The average value of f>j_,BBL> c /- Fig. 1131 over im- 
pact parameters of the proton-proton collisions, Eq. 15711 . as 
a function of the resolving parton's momentum fraction, xr. 
Shown are the averages computed with the impact parame- 
ter distribution corresponding to the hard dijet trigger, P2 
(solid line), and the double dijet trigger, P4 (dashed line), 
for yfs = 14000 GeV and dijet momentum p± = 25 GeV, cf. 

Fig. n 



partons is well justified. What is equally important, the 
6-distributions implied by the hard multijet trigger sup- 
press the contributions from large impact parameters, 
b > 1 fm, where p\ BBL drops below ~ 1 GeV 2 , mean- 
ing that the gluon density seen by the small-x partons 
is so low that the BBL is never reached in the region 
where our approximations are justified. For the above 
dijet trigger, the fraction of events with b so large that 
P± bbl(^) < 1 GeV 2 is no larger that — 10% for x R > 0.1; 
for the double dijet trigger it drops to ~ 1%. In contrast, 
for generic inelastic events without the hard multijet trig- 
ger, the 6-distributions of Fig. show that a significant 
fraction of events would involve 6-values for which our 
perturbative treatment would not be justified. In this 
sense, the possibility to control the effective impact pa- 
rameters in the pp collisions by way of the hard multijet 
trigger is more than just an enhancement of otherwise 
well-defined contributions — it is crucial for the very 
applicability of the estimates presented here. 

The above estimates are based on the interaction of 
color-octet (88) dipoles, corresponding to gluon partons, 
with the gluon field of the other proton. They can eas- 
ily be extended to the case of color-triplet (33) dipoles, 
corresponding to quark partons, see Eq. and after. 
In particular, one finds that in this case the values of 
(pj_ BBL ) are approximately 0.5 times the values for glu- 
ons over the range of xr shown in Fig. 1141 for both the 
dijet and the double dijet trigger. 



It is worth noting that our estimate allows to avoid 
double counting due to multiple rescattering at the vir- 
tualities much smaller than the BBL. At the same time, 
it neglects an additional broadening due to the hard scat- 
tering at scales somewhat larger than BBL. 



VI. FINAL STATE PROPERTIES FOR 
CENTRAL pp COLLISIONS 

We have seen that the probability for partons in one 
proton to interact with the other proton near the BBL is 
greatly increased by a trigger on hard dijet production, 
which acts as a centrality filter. We now outline the 
consequences of such interactions close to the BBL for 
the final state of the hadronic collision. 

Generally speaking, partons propagating through a 
strong gluon field get transverse momenta of the order 
of the maximum transverse momentum for which the 
interaction remains black, p_l,bbl cf. the discussion in 
Refs. |2(| I27I l28l |29| . This value of transverse momen- 
tum is (approximately) related to the maximum Q 2 at 
which the BBL is still valid as 



P_l,bbl ~ Qbbl/2, 



(58) 



assuming that in the hard interaction a system with mass 
of the order of Q is produced. This relation can also be 
expressed as 



P_L,BBL 



2d 



3 

2d' 



(59) 



where d is the dipole size, cf. Eq. (|47|l . It is interesting 
that a similar numerical estimate follows from the uncer- 
tainty relation: Regarding d as conjugate to p± of the 
parton in the dipole one obtains p^.bbl ~ 7r/(2d). Note 
that our definition of the BBL scale, Qbbl, is different 
from the saturation scale Q s of the Color Glass Conden- 
sate picture; nevertheless the two scales are numerically 
comparable. 

We have seen that in the central collisions selected 
by the "centrality trigger" the maximum transverse mo- 
menta squared, p\ BBL , are 3> lGeV 2 , cf. Fig. El This 
allows us to describe the hadronization of these partons 
assuming independent fragmentation. In this approxima- 
tion, the differential cross section for the semi-inclusive 
production of a given hadron h, characterized by its ra- 
pidity relative to the parent proton of the resolving par- 
ton, logz, and transverse momentum , fcj_, is given by 

m 



da pp - 



dz d 2 k± 

- E 

x 6^(1^ 



d 2 q± d 2 l ± 



dx — 
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X 



fi/pi^i 1 



[z/x)qj 
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>J ^ D i/p (z/x,l±,Q 2 



BBL 



)<q±). (60) 
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Here fi/ p (x,Q 2 ) is the leading-twist parton (i = quark, 
gluon) distribution in the proton, c(q±) the transverse 
momentum distribution of the partons after passing 
through the strong gluon field, which can be estimated 
in various models (see e.g. Ref. H3|), and D i j h the par- 
ton fragmentation function. Here we make the natural 
assumption, in line with the discussion above, that the 
small-x partons are resolved at the scale Q = Q| BL , 
i.e. that the factorization scale is given by the maxi- 
mal virtuality at which the interaction is close to the 
BBL. In Ref. [2^ a similar approach was applied to par- 
ticle production in central proton — nucleus collisions. We 
stress that the approximation of independent fragmenta- 
tion is justified for sufficiently large values of z (leading 
hadrons); the relevant range in z grows with the trans- 
verse momentum of the produced hadrons. For small 
z coherence effects need to be taken into account, see 
Ref. for a discussion. 

In the BBL, estimates based on Eq. ljl)D)) should be 
regarded as an upper bound for the spectrum of lead- 
ing particles, since in the case of large-angle scattering 
a parton may actually convert into two high-p^ partons 
with, on average, equal light cone fractions, resulting in 
an even steeper drop of the leading hadron spectrum with 
z. This is what happens in deep -inelastic scattering of 
a virtual photon in the BBL [26j . In our case we expect 
somewhat smaller suppression due to this effect, because 
large impact parameters contribute to the process (al- 
though with small probability), and in this case much 
more forward particles would be produced. At the same 
time, the abovementioned trend of a steeper drop of the 
spectrum with z should be much more pronounced for 
large-fcj_ particles, since in collisions at large impact pa- 
rameters very few large-A:^ particles would be produced. 

A detailed numerical investigation of hadron produc- 
tion in pp collisions based on Eq. (|tJ0|> is beyond the scope 
of the present paper. Here we only list some expecta- 
tions for the qualitative properties of hadron production 
at large rapidities which follow from independent frag- 
mentation in the BBL: 

• The leading particle spectrum will be strongly sup- 
pressed compared to interactions far from the BBL. 
The suppression will be especially pronounced for 
nucleons, so that for z > 0.1 the differential multi- 
plicity of pions should exceed that of nucleons. 

• The average transverse momenta of the leading par- 
ticles will be > 1 GeV/c. 

• There will typically be no correlation between the 
transverse momenta of leading hadrons, since they 
originate from two different partons which have un- 
corrected transverse momenta. Some correlations 
will remain, however, because two partons pro- 
duced in collisions of small-a; and large-a; partons 
may end up at similar rapidities. 

• For small impact parameters (which constitute a 
relatively small fraction of the total inelastic cross 



section but dominate in new particle production) 
a large fraction of the events will have no particles 
with z > 0.02 0.05. This suppression will occur 
simultaneously in both fragmentation regions, cor- 
responding to the emergence of long-range rapid- 
ity correlations between the fragmentation regions. 
For studies of this feature of the central pp collisions 
it would be desirable to have good acceptance for 
both leading charged and neutral particles. This 
would allow one to measure the fraction of events 
without leading particles as a function of the cen- 
trality of the collision. 

• In the forward production of dimuons or dijets one 
expects a broadening of the distribution over trans- 
verse momenta |3l|. as well as a weaker depen- 
dence of the dimuon production cross section on 
the dimuon mass for masses < few GeV 29]. 

• The background for heavy particle or high-pj_ jet 
production should contain a significant fraction of 
hadrons with transverse momenta p±_ ~ Pj_,bbLj 
originating from fragmentation of partons affected 
by the strong gluon field. The direction of the 
transverse momenta of these hadrons should be un- 
related to the transverse momenta of the jets. This 
phenomenon will make it difficult to establish the 
direction of jets unless p± (jet) ^> p_l,bbl- 

In connection with particle production at large rapidities, 
another interesting quantity to study will be the incident 
energy dependence of the leading particle xf multiplicity 
with the dijet "centrality trigger" for fixed values of x\ 
and X2- In the absence of the trigger this would corre- 
spond to the study of usual Feynman scaling violations. 
In this case it is known that different impact parame- 
ters contribute, and that leading particle production is 
predominantly a large impact parameter phenomenon, 
which likely leads to a rather weak violation of Feyman 
scaling. On the other hand, in the "conditional multiplic- 
ity" with the centrality trigger small impact parameters 
give the dominant contribution, and thus the suppression 
of the forward spectrum should strongly increase with en- 
ergy. 

The proximity to the BBL in central pp collisions will 
also lead to observable effects in particle production at 
small rapidities. One expects a significant increase of 
the multiplicity at small rapidities, because interactions 
in the BBL will likely lead to generation of large color 
charges in the fragmentation regions. Also, as we dis- 
cussed in Section llVI the production of multiple minijets 
will be strongly enhanced. Such an increase should in 
fact be present already at the Tevatron collider, in events 
with a trigger on two-jet or Z° production. We are aware 
of only one study [3^] which investigated the correlation 
of the underlying event structure with the presence of 
such a trigger. An increase of the multiplicity at small 
rapidities was indeed observed. 

The detailed modeling of the phenomena outlined in 
this section will require building a Monte Carlo event 
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generator, accounting for the proximity of the spectator 
interactions to the BBL and for the new pattern of flow 
of color excitations. 

VII. CONCLUSIONS 

In this paper we have demonstrated that a trigger on 
hard dijet production strongly reduces the effective im- 
pact parameters in pp collisions at LHC and, to some 
extent, at Tevatron energies. The possibility to select 
central collisions with a reasonable rate is of consider- 
able practical importance. 

We have argued that the structure of the final states 
in central pp collisions at LHC energies will differ sig- 
nificantly from that of minimal bias events. The reason 
is that in the central transverse region the interaction 
of hard partons with the gluon field in the other pro- 
ton approaches the unitarity (black-body) limit, leading 
to an enhancement of transverse momenta and depletion 
of longitudinal momenta in hadron production at large 
rapidities. 

The proposed centrality trigger offers new opportuni- 
ties for realistic studies of the physics of strong gluon 
fields in QCD. On the experimental side, it singles out 
central events for which the chances of reaching the BBL 
are maximized, and its signatures in the final state can 
be clearly identified. On the theoretical side, it quan- 
titatively defines the region in transverse space in which 
one should expect deviations from DGLAP evolution due 
to unitarity effects. This information should be incorpo- 
rated in studies of non-linear QCD evolution within the 
renormalization group approach, see e.g. Ref. |3^ | and 
references therein. So far, such studies have assumed 
infinite extension of the "dense gluon medium" in the 
transverse plane. They have also neglected the contribu- 
tions of large-x partons to the parton densities at small 
x due to log Q 2 -evolution, which are taken into account 
in our approach. 

New heavy particles at LHC will be produced practi- 
cally only in central collisions, as can be selected with the 
proposed hard dijet trigger. The fact that such collisions 
are also the ones in which the interactions of hard partons 
reaches the BBL poses new challenges for the analysis of 
such events. One needs to identify the signatures of new 
particle production on top of the very specific modifica- 
tions of the final state implied by the BBL in the central 
region. This could imply e.g. changes in the definition 
of jets due to the enhanced pedestal of soft hadrons and 
enhanced production of minijets, as well as changes in 
the cuts necessary to define an isolated lepton. Further- 
more, the sizable "intrinsic" p± of the partons acquired 
from interactions with the strong gluon field may impose 



limits on the accuracy of the determination of the masses 
of produced heavy particles. 

There is an interesting connection of the pp collisions at 
LHC energies discussed here with cosmic ray physics near 
the Greisen-Zatsepin-Kuzmin (GZK) cutoff, see Ref. [34| 
for a recent review. The density of gluons through which 
a proton propagates in a central pp collisions at LHC 
energies is comparable to the typical density encountered 
by a proton near the GZK cutoff in collisions with air 
(i.e., light nuclei). Thus, the BBL effects on forward 
particle production described above could have an impact 
on the energy spectrum and composition of cosmic rays 
near the GZK cutoff. 

Learning how hadron production depends on the im- 
pact parameter of the pp collision will allow to address 
also other questions of strong interaction dynamics, not 
primarily related to the BBL. For example, it is often 
argued that events with large hadron multiplicities (a 
factor of > 2 larger than average) at small rapidities are 
due to collisions at small impact parameters, in which 
the soft parton clouds of the colliding protons overlap 
much stronger than in average collisions. Once the cor- 
relation between the centrality of the collision and the 
suppression of the forward spectra etc. has been estab- 
lished with the help of the dijet trigger, one would be 
able to check whether similar effects are present also in 
the high-multiplicity events, testing the hypothesis that 
these are central collisions. 

Finally, the comparison of our model estimate of the 
cross section for double dijet production with the CDF 
data indicates that there may be significant spatial 
correlations of partons in the transverse plane. This 
interesting phenomenon should be investigated further. 
Not only is it important for possible extensions of the 
single dijet trigger to multiple dijet production — it may 
also reveal interesting information about the structure 
of the nucleon at low scales, which generates the parton 
distributions at the dijet production scale through Q 2 - 
evolution. 
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